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We apply density functional theory, in the local density approximation, to a quasi-one-dimensional 
electron gas in order to quantify the effect of Coulomb and correlation effects in modulating, and 
therefore patterning, the charge density distribution. Our calculations are presented specifically for 
surface-gate-defined quasi-one-dimensional quantum wires in a GaAs-AlGaAs heterostructure but 
we expect our results to apply more generally for other low dimensional semiconductor systems. We 
show that at high densities with strong confinement, screening of electrons in the direction transverse 
to the wire is efficient and density modulations are not visible. In the low-density, weak-confinement 
regime, the exchange-correlation potential induces small density modulations as the electrons are 
depleted from the wire. At the weakest confinements and lowest densities, the electron density 
splits into two rows thereby forming a pair of quantum wires that lie beneath the surface gates. An 
additional double-well external potential forms at very low density which enhances this row splitting 
phenomenon. We produce phase diagrams that show a transition between the presence of a single 
quantum wire in a split-gate structure and two quantum wires. We suggest that this phenomenon can 
be used to pattern and modulate the electron density in low-dimensional structures with particular 
application to systems where a proximity effect from a surface gate would be valuable. 


I. INTRODUCTION 

Very large scale integration and the invention of the 
microprocessor has revolutionized almost every aspect 
of the modern world. The importance of the miniatur¬ 
ization of the field effect transistor (FET) is most dra¬ 
matically demonstrated by the celebrated Moore’s Law, 
whereby the size and cost of transistors decreases at an 
exponential rate. However, since its inception, it has 
been known that Moore’s Law will inevitably fail at the 
atomic level as quantum effects become dominant and the 
conventional operation mode of FETs breaks down. Un¬ 
derstanding and manipulating quantum effects in these 
devices will play an essential role in extending Moore’s 
law into the twenty-first century. 

In semiconductor heterostructures, the FET is realised 
through a combination of band-gap engineering and elec¬ 
trostatic gating. Dopants are used to promote elec¬ 
trons into the material’s conduction band whilst pat¬ 
terned metallic surface gates define electrostatic poten¬ 
tials which alter the density and the dimensionality of the 
electron system niia. Using these techniques, the elec¬ 
trons can be confined to structures with features which 
are on the same length scale as the surface gate features. 

In conventional devices, engineered electrostatic po¬ 
tentials are the only method for confining the electrons. 
However, at low densities electron screening becomes less 
efficient and Coulombic repulsion can be used in con¬ 
cert with electrostatic confinement to generate features at 
smaller length scales. For example, whilst single electrons 
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can be localised into a set of coupled quantum dots using 
sophisticated gate arrays EHS], spontaneous localisation 
of the electron gas into a Wigner crystal state provides a 
simpler alternative for generating an ordered set of elec¬ 
trons. Quantum Monte Carlo [6] and density functional 
theory cm simulations have shown that, within a quasi- 
one-dimensional system, interactions cause the electrons 
to localise and form a one-dimensional Wigner crystal. 
As the confining potential is weakened, the crystal under¬ 
goes a phase transition first to zig-zag, then double row 
configurations [MI]. Experimental evidence of the tran¬ 
sition to a localised electron state has been observed in 
electron transport measurements in GaAs/AlGaAs het¬ 
erostructures. A quasi-one-dimensional wire can be de¬ 
fined using metallic surface gates iiaiis] and an addi¬ 
tional top gate can be used to vary the electron density in 
the wire pAUTB] . In the high density regime, the electrons 
behave as non-interacting particles and the conductance 
is quantised in units of 2e^/h [17] but as the electron’s 
density in the wire is reduced, the first plateau suddenly 
disappears. Further reduction in the density leads to a 
revival of this plateau. It has been speculated that this 
“missing plateau” can been attributed to the degeneracy 
of the zig-zag state m- The electron separation in these 
devices is considerably smaller than could be achieved 
using surface gate patterning. 

In this paper, we propose a method for exploiting 
Coulombic repulsion in order to create two quantum 
wires in the electron gas with a separation width sim¬ 
ilar to the length scale of the surface gate pattern. We 
use density functional theory (DFT) in the local density 
approximation (LDA) to find the ground-state electronic 
structure of translationally invariant quantum wires with 
experimentally realistic, transverse confining potentials. 
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The confining potentials are calculated exactly, includ¬ 
ing all anharmonic and image charge contributions, by 
solving Poisson’s equation using finite-element analysis 
(FEA) with adaptive gridding [18] . These computational 
techniques allow us to probe a regime where electrons 
are confined through a combination of electrostatic and 
Coulombic potentials, which opens the door to a wide va¬ 
riety of applications in future quantum one-dimensional 
technologies. 

We start by showing that in the strong-confinement 
regime, electron screening almost completely removes 
transverse density modulations. These results contradict 
the observations of Laux et al. [2] and are due to the 
highly accurate FEA potentials used. We extend these 
results into the low-density, weak-confinement regime 
and find that transverse density modulations exist at low 
densities for a range of confinements. These modulations 
are only observed in the EDA and therefore are driven 
by exchange-correlation effects. For the weakest possi¬ 
ble confinement and the lowest achievable densities, the 
electrostatic repulsion of the electrons splits the density 
into two rows. For the experimental setup used to reach 
this regime, it is possible for an external double-well po¬ 
tential to form as the top gate is negatively biased with 
respect to the side gates. Our simulations show that this 
does occur but only after a double-row configuration of 
the density has appeared due to Coulombic interactions. 


II. MODEL 

Semiconductor heterostructures consist of layers of 
semiconductor material with differing alloy composition, 
typically grown by molecular beam epitaxy, which pro¬ 
duce a variable band-gap material. Ionized dopants in¬ 
troduced into specific layers generate internal electric 
fields which can bend the conduction band at a remote 
interface below the chemical potential and allow low¬ 
dimensional electron systems to form. Additional fields 
can be applied using metallic surface gates to shape the 
electrons into quasi-one-dimensional quantum wires. In 
order to accurately calculate the structure of the electron 
density for these low-dimensional systems, we will include 
the effects of dopant densities, surface charges and biased 
metallic surface gates on arbitrary layered semiconduc¬ 
tor heterostructures. By simulating the potential from 
the entire device, and not resorting to analytic models, 
anharmonic contributions present in the confining poten¬ 
tial will be included which can strongly affect the form 
of the electron density. 

Firstly, we must calculate the electrostatic potential 
0(r). For a given charge density p(r), (/)(r) is the solution 
to Poisson’s equation 

g p(rl] ^ -V • (e(rlV(A(rl) -^=0 (1) 

^0 


As the rows decouple, the energy difference between 
the ground and first excited sub-bands of the wire van¬ 
ishes and the first conductance plateau disappears. It is 
possible that disorder could break this degeneracy as the 
confinement potential weakens, which would lead to the 
revival of this plateau as measured in Refs. [MIITG] . but 
this effect is not included. The results demonstrate that 
anharmonic contributions to the confinement potential 
are crucial to understanding the electronic structure of 
weakly-confined one-dimensional wires. Electron locali¬ 
sation in this regime is strongly determined by the inter¬ 
play between electron-electron and external confinement 
potentials. 


The paper is organized as follows: In Sec. [11] we de¬ 
scribe the model used to calculate the ground states den¬ 
sity of quasi-one-dimensional quantum wires with arbi¬ 
trary semiconductor heterostructures and surface gate 
geometries. We use this model in Sec. Ill to show that 
in the strong-confinement regime electron screening is ef¬ 
ficient and density modulations effectively disappear. In 
Sec. |IV[ we explore the electronic structure in the low- 
density, weak-confinement regime and we report the vari¬ 
ous phases of the quasi-one-dimensional wire. We present 
our conclusions in Sec. [Vj along with possible experimen¬ 
tal methods for measuring the density modulations pre¬ 
dicted in this paper. 


where e(r) is the relative permittivity [19] of the layered 
semiconductor materials at r = (x^y^z). The bound¬ 
ary conditions on 0(r) are determined by the voltages on 
the metallic surface gates. For free surfaces, we apply 
von Neumann boundary conditions n • V0(r) = 0 where 
n is the normal to the domain boundary. We calculate 
the potential using finite-element analysis (FEA) on an 
adaptively refined mesh using the deal.II finite-element 
library [18] which increases the computational efficiency 
and accuracy of the solution as well as eliminating high- 
order Fourier components in 0(r) which may be gener¬ 
ated by the discretisation onto a regular grid. 

Secondly, we need to calculate the charge density p{f) 
for an arbitrary semiconductor heterostructure with an 
electrostatic potential 0(r). We start by separating p{f) 
into static and mobile charges. Our method for calcu¬ 
lating the static charges - which consist of the dopant 
and surface charges - is described in App.[^ The mobile 
charges considered in this paper are electrons confined to 
a narrow region below a GaAs/AlGaAs interface. To cal¬ 
culate this density, we need to solve the quantum many- 
body problem, which is impossible, so approximations 
need to be made. The approach used in this paper is to 
use density functional theory (DFT) within the Hartree 
and local density approximations (LDA). 

The electron charge density pe{r) = —en{f) = 
is calculated by solving the Kohn- 
Sham (KS) equations in the effective mass approxima- 
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tion 


Stern [21], we assume that 


2m^ 


+ K.t(r) + VH[n]{f^ + 14c[n](r-)j 


where = O.OGlrrie is the reduced mass of the electron 
in GaAs. The external potential Vext{r) generated by the 
dopants, surface charges and applied metallic gate bias 
voltages and the classical electron-electron or Hartree po¬ 
tential VH[n]{f) are both obtained from the solution to 
Eq.[T] 


Vxt(r) + Vh[ n ](rl = (rl - e^{f^ (3) 


dpe{(f>nii^) 

d(t>n{i”) 


j 

d 

■—J f{E,T)g^D{E-,ej)dE 


( 7 ) 


The Newton step is solved using the same TEA methods 
as Eq. The mixing parameter tn is calculated by find¬ 
ing the zero of the function F{t) = {(t)n + tx^) • x^ 

which accelerates convergence [22]. We update the 
exchange-correlation potential using a simple function 
mixing method 


14 c,newN(f) = (1 - ^)E,e,old[n](f) + ^^^c,calc N(r) (8) 


where, for layered semiconductor heterostructures, the 
band gap Eg{r) is a function of the growth direction 
^ only. Vxc[^]{f) is the exchange-correlation potential 
which contains all quantum many-body effects. Eor the 
Hartree approximation Vxc[^]{r) = 0 whilst in the EDA 
we use the parameterization of Perdew and Wang [20] in 
effective atomic units with 1 Ry* = and 1 

~ h‘^eoer/m^e‘^. 

One-dimensional devices are translationally invariant 
in the longitudinal direction x so the solutions to the 
Kohn-Sham equation are plane waves 

(4) 


where a = 0.3 is a constant mixing parameter. The 
iteration scheme terminates when the maximum change 
in the potential V{y^z) for the grid on which the electron 
density is calculated is less than 0.1/reV. 

To ensure that the converged values of (j){f) and n(r) 
are consistent, the density obtained using this iterative 
procedure is used to calculate a new initial potential 0o- 
Only when the maximum value of the initial Newton step 
xq is less than 4/ieV do we consider the density to have 
reached convergence. 

III. ELECTRON DENSITIES IN THE STRONG 
CONFINEMENT REGIME 


Inputting Tj (r) into Eq. gives 


+ +'E(y,d )V’i(2/,d 


2m* \dy‘^ dz^ 


= Ej'ipj{y,z) (5) 


where V{y,z) = Eg{z)f2 - e(j){y,z) + Vxc[n]{y, z). The 
transverse wave functions 'iljj{y^z) are independent of k 
so the parabolic dispersion relation for each of the one¬ 
dimensional sub-bands j can be integrated numerically 
to give the electron density: 


/ C 

-c 


f{E,T)g,n{E;Sg)\ijgiy,z)fdE (6) 


where ej = Ej — k?k^/2m^^ — 

^/2m^l{E — Sj) is the one-dimensional density of 
states and f{E,T) = (1 + ^ElksT^-i Eermi-Dirac 

function for the system at temperature T. 

Einally, we use an iterative procedure to find a con¬ 
sistent solution for Eqs.[^and[^ An initial electrostatic 
potential 0o is calculated for the system with n(r) = 0. 
This potential is updated using an adaptive Newton 
step method 0n+i = where tn is a mixing 

parameter and the Newton step Xn is the solution to 
9'{4>n)xn = where g{(pn) = 9 [(f>n{r), p{(pnX)){r)] 

and g'{(l)n) is the Jacobian of g{(l)n)- Eollowing Laux and 


Eor a truly one dimensional system, electrons are con¬ 
fined so that only motion in the longitudinal direction 
is possible. However, in any realistic physical system, 
the confining potential is finite and the electrons will be 
able to move in the transverse direction which defines a 
quasi-one-dimensional wire. We can define a strong con¬ 
finement regime, where the electrons are forced into the 
center of the wire and there is no irregular structure of the 
electrons. The electrons will spread out to some extent 
in the transverse direction as the density increases but, 
when the confinement is strong, screening is expected to 
be efficient so any density modulations should be washed 
out. 

In semiconductor heterostructure devices, a pair of 
parallel negatively-biased metallic gates above an intrin¬ 
sic two-dimensional electron gas will deplete the elec¬ 
trons underneath and confine them into a quasi-one- 
dimensional quantum wire. Such a device was modelled 
by Laux et al. [2] and is described in detail in App. 

In this section, we reproduce their results and show that 
screening of the electron density is efficient in these de¬ 
vices. 

Cross-sections 6nm below the GaAs/AlGaAs interface 
of the self-consistent potential V{y, z) for Hartree and 
LDA calculations are shown in Eig. [^ When the wire 
is occupied, the potential is essentially flat, which shows 
that electron screening is efficient in this device. The elec¬ 
tron density along this same line is shown in Eig. [^ As 
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Lateral Distance (nm) Lateral Distance (nm) 

(a) Hartree (b) LDA 

FIG. 1. (Color online) Lateral potential at 4.2K on a cross-section 6nm below the GaAs/AlGaAs interface for the device 
described in Sec. IIIII 



(a) Hartree 



(b) LDA 


FIG. 2. (Color online) Electron charge density at 4.2K along the same cross-section as in Fig. The electron screening is 
efficient and density modulations in the transverse direction vanish for this strong conhnement regime. The difference between 
the Hartree and local-density approximations is minimal. 


observed in the original paper, the exchange-correlation 
potential does not have a significant effect. Laux et al 
saw oscillations in the electron density but these do not 
appear in our calculations. Although we were not able 
to reproduce their modulations, discretizing the poten¬ 
tial on a regular grid could have prevented screening at 
length scales similar to the lattice spacing. We believe 
that the disappearance of the density modulations is due 
to the high accuracy potentials calculated using the FEA 
method. 


IV. ELECTRON DENSITIES IN THE 
LOW-DENSITY, WEAK-CONFINEMENT 
REGIME 


We now turn our attention to the low-density, weak- 
confinement limit of the quasi-one-dimensional wire. The 
split gate geometry from Sec. [ITT] does not allow the den¬ 
sity and confinement of the wire to be tuned indepen¬ 
dently. As the split gate voltage decreases, the transverse 
confining potential Wxt(^) at the GaAs/AlGaAs interface 
in the ^-direction becomes weaker but its minimum de¬ 
creases. Electrons can flow into the system and the den¬ 
sity increases. The top-gate, split-gate device used in the 
series of experiments initiated by Hew et al. [HHin] al- 
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lows the confinement and density to be tuned separately. 
The surface of the split gate device is covered with a 
layer of insulating polymer on top of which an additional 
metallic gate is deposited. The voltage on the top gate 
Vtg can be varied separately from the split gate voltage 
Vsg such that the confining potential 14xt can be weak¬ 
ened whilst keeping the electron density constant. In this 
section, we look at such a device which is described in 
more detail in App. |B 2[ The electron density was calcu¬ 
lated using the effective-band approximation described in 
App.[^which produces very good approximate solutions 
to the Kohn-Sham equations. 

Self-consistent solutions for this system were calculated 
for split gate voltages Vgg in the range -0.58V to -0.46V 
and top gate voltages Vtg in the range -1.75V to -I.IOV. 
This range of Vgg corresponds to the weakest possible 
confining potential as the electrons are no longer depleted 
under the split gate for Vsg > —0.46V. The range of Vtg 
corresponds to the low density regime where the con¬ 
striction is almost entirely drained of electrons. It is 
in this regime that we expect electrostatic interactions 
and exchange-correlation effects to be enhanced. In this 
section, we will describe the electronic structure of the 
quasi-one-dimensional quantum wire in this regime. 

For a fixed split gate voltage, as the top gate volt¬ 
age becomes more negative and the density decreases, 
there is a qualitative difference induced by the exchange- 
correlation potential. In the LDA calculations, we ob¬ 
served density modulations in the ^-direction transverse 
to the wire axis which do not appear in the Hartree 
approximation, where only electrostatic interactions are 
taken into account. A comparison of the different be¬ 
haviours for Hartree and LDA calculations is shown in 
Fig. 1^ for Vsg = —0.50V. The varying top gate volt¬ 
age changes the slope of the conductance band so the 
maximum of the electron density in the growth direc¬ 
tion shifts. Therefore, for the top-gate, split-gate devices 
simulated in this section, the density has been integrated 
over the ^-direction. For large densities, the electrons 
effective screen themselves but in the LDA calculations, 
as they are depleted from the wire, the modulations in 
the transverse density distribution appear. The den¬ 
sity modulations observed here are qualitatively different 
from those observed in the papers by Laux et al [2] and 
Malet et al. [23] as they are produced exclusively by the 
exchange-correlation potential. 

In the weak-confinement regime, for high densities the 
electrons also screen themselves effectively but as the 
electrons are depleted, the charge distribution localises 
to the edges of the wire and a double-row density distri¬ 
bution forms. Due to the translational invariance of our 
simulations, this density distribution corresponds to both 
the zig-zag and double row phases of the weakly confined 
quantum wire m- For most of the voltage values Vsg and 
Vtg where there are electrons in the wire, the external po¬ 
tential Vext{r) at the GaAs/AlGaAs interface has a sin¬ 
gle minimum in the transverse direction. Row formation 
starts to occur in this regime where it is driven purely by 


the electrostatic repulsion of the electrons. However, for 
the most negative values of Vtg^ the top gate generates a 
double-well potential with two minima in the transverse 
direction. This generates an additional anharmonic force 
on the electrons which is independent of the many-body 
effect. 

A typical set of integrated transverse densities in this 
weak-confinement regime is shown in Fig. for Vsg = 
—0.46V. Whilst there is some additional localisation for 
the LDA calculations, the main feature is the splitting of 
the electron density into two rows. As the electron den¬ 
sity is reduced, this is driven first by the electrostatic re¬ 
pulsion of the electrons and then by the external double¬ 
well potential. Before the electrons are completely de¬ 
pleted from the wire, the two rows are practically decou¬ 
pled and the device consists of a pair of quantum wires 
which lie beneath the surface gates. 

The number of peaks in the transverse density distribu¬ 
tion is shown in Fig. as a function of Vtg and Vsg. The 
boundary where Wxtv^) just below the GaAs/AlGaAs 
interface becomes a double-well potential is denoted by 
the green dotted line. We resolve the peaks in the den¬ 
sity using a valley definition; two peaks exist when the 
minimum value of the density between the two peaks 
is less that a given factor of the maximum density, ie. 
'^min < This definition allows us to exclude spu¬ 

rious peaks when the density is approximately constant 
across the width of the wire. For this paper, we have 
chosen /3 = 0.99. The double-row regime is well defined, 
with a boundary which is similar for both Hartree and 
LDA calculations. The transverse density modulations 
are observed for stronger confining potentials just before 
the quantum wire is depleted of electrons and only occur 
in the LDA calculations where the exchange-correlation 
potential is non-zero. For large densities, the electrons 
screen themselves efficiently and there is no fine struc¬ 
ture in the electron density distribution. 

The voltages Vsg and Vtg determine the density of elec¬ 
trons and the confinement in the transverse direction 
but there is no simple relation between these voltages 
and the physical parameters of the quasi-one-dimensional 
electron system. We can express the data in Fig. in 
terms of the two physical length scales of the system: 
the Wigner-Seitz radius r^, defined here as the aver¬ 
age separation of the electrons in the longitudinal di¬ 
rection, and the distance in the transverse direction at 
which the electrostatic and confinement energies are the 
same tq. For a model harmonic potential of frequency 
cj, ro = (2e^/eoe^m*ci;)3 [T0l[24]. However, in our sim¬ 
ulations, the electrostatic potential Vext{'^) contains an¬ 
harmonic contributions. Therefore, we define ro as the 
distance for which 

Vextiro) = e^/47rereoro (9) 

For a double-well potential, ro is defined as the separation 
between the minima of the potential if this is greater than 
the value of ro calculated from Eq. The lengths r^ and 
ro are normalised by the effective Bohr radius for GaAs 
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Lateral Distance (nm) Lateral Distance (nm) 

(a) Hartree (b) LDA 

FIG. 3. (Color online) Integrated transverse density for Vsg = —0.50V at 50mK. The electrons do not form a double-row 
density but the exchange-correlation potential generates small transverse density modulations which do not occur in the 
Hartree simulations. 




Lateral Distance (nm) Lateral Distance (nm) 


(a) Hartree 


(b) LDA 


FIG. 4. (Color online) Integrated transverse density for Vsg = —0.46V at 50mK. Traces show the density for different top-gate 
voltages from Vtg = —1.70V to Vtg = —1.40V in steps of AVtg = 0.02V. This is the weak-conhnement regime where electrostatic 
repulsion is dominant and double-row formation occurs. For high densities, the electrons screen themselves but as the density 
decreases, two rows forms in both the Hartree and LDA calculations. For AVtg < —1.58V, the external potential Vext(^) 
acquires a double-well structure which enhances row separation. When the wires are completely decoupled, the electron density 
lies almost entirely underneath the surface gates. 


The number of peaks in the transverse density distri¬ 
bution is plotted in Fig as a function of Vg and ro for 
the LDA calculations. We use the Hartree calculations 
to define the boundary between the exchange-correlation 
potential driven density modulations and double-row for¬ 
mation, which is denoted by the blue dotted line. For 
small ro, there is a strong confining potential and the 
density modulations are screened out in the Hartree ap¬ 
proximation. However, for the LDA, as the confining 
potential is weakened, it is possible to generate a multi- 


peaked density at low density. For large ro, the confining 
potential is weak and electrostatic repulsion forces the 
electrons into a double-peaked density distribution. Such 
a density is consistent with a correlated system where the 
electrons are arranged in either a zig-zag configuration 
or in two parallel rows. Interestingly, for lower densities 
the electrostatic transition to the double-row density oc¬ 
curs at weaker confining potentials, which is consistent 
with the quantum Monte Carlo calculations of Mehta et 

al. [To]. 

The number of Kohn-Sham bands below the chemical 
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FIG. 5. (Color online) Number of peaks in the transverse density as a function of Vtg and Vsg. For weak confinement and low 
density, double-row formation occurs for small Vsg and large Vtg in both the Hartree and LDA calculations. Transverse density 
modulations in the LDA calculations at low density create a peaked structure which extends to a less weakly conhned regime 
where Vsg is more negative. Below the dashed green line, the I4xt(^) forms a double-well potential below the GaAs/AlGaAs 
interface. 



0 5 10 15 20 

ro(nm) 

FIG. 6. (Color online) Number of peaks in the transverse 
density as a function of ro and Vs in the local density approx¬ 
imation. A typical split gate devices, for example in Sec. |IIH 
corresponds to the small Vs and ro limit. The blue dotted line 
dehnes the transition above which the electrons separate into 
two rows in the Hartree approximation which corresponds to 
zig-zag or double row conhgurations and is due to electrostatic 
repulsion. Modulations in the transverse density generated 
by the exchange-correlation potential occur for larger Vs and 
ro. Calculations where Ve^t{r) has two minima correspond to 
ro > 10 and are not shown in this hgure. 


potential is shown in Fig. In the Landauer-Biittiker 
formalism [Tans], this is equivalent to the conductance 
of the quasi-one-dimensional wire G in units of 2e^/h. 
The plateaux edges are infinitely sharp as the wire is 
translationally invariant in the transport direction. In 


both the Hartree and LDA approximations, we see that 
for a fixed top gate voltage Vtg > —1.6V, the conduc¬ 
tance is quantised in units of 2e^/h until the wire def¬ 
inition is lost at Vsg = —0.46V. For more negative top 
gate voltages, we see a jump from zero to two occupied 
sub-bands. As the electron rows separate, the coupling 
between the symmetric and antisymmetric Kohn-Sham 
bands tends to zero. This degeneracy leads to an initial 
jump in the conductance to 4e^//i as the rows are de¬ 
coupled and form two independent quantum wires. In 
our simulations, there is no mechanism for breaking this 
degeneracy and regaining the missing plateau at 2e^ jh 
as seen in Refs. pMHT6] . A possible explanation for this 
effect would be a disorder-driven breaking of the parity 
symmetry of the system in the ^-direction. However, the 
inclusion of disorder, from the ionized dopant background 
for example, is outside the scope of this paper. 

V. CONCLUSIONS 

In this paper, we have studied the ground-state elec¬ 
tronic structure of the quasi-one-dimensional quantum 
wire. Self-consistent electron densities and electrostatic 
potentials for realistic semiconductor heterostructures 
were calculated with density functional theory using 
finite-element analysis. We looked at the strong confine¬ 
ment regime by simulating a simple split-gate device, re¬ 
producing the simulations of Laux et al [2] , using FEA to 
calculate accurate electrostatic potentials. We saw that 
the transverse density modulations previously observed 
were screened out by the electron-electron interactions. 

We extended these results into the weak-confinement, 
low-density regime by simulating a split-gate, top-gate 
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Vsg{V) Vsg{V) 

(a) Hartree (b) LDA 

FIG. 7. (Color online) Number of Kohn-Sham bands below the chemical potential as a function of Vtg and Vsg. A double-jump 
in the number of conducting bands is seen for very weak-confinement and is driven by the formation of a double-well potential 
in the external gate potential I4xt(^) below the dashed green line. The conduction bands do not correlate with either the 
formation of a double-row in the density or the density modulations generated by the exchange-correlation potential. 


device, which allowed us to control the density and trans¬ 
verse confinement of the quantum wire separately. We 
showed that exchange-correlation effects produce small 
modulations in the transverse electron density for low 
electron densities across a range of confining potentials. 
In the weak-confinement, low-density limit, the electrons 
separate into two rows due to electrostatic repulsion. For 
the realistic experimental device we simulated, an exter¬ 
nal double-well potential forms which enhances this sep¬ 
aration after row formation has occurred. Decoupling of 
the rows leads to a degeneracy in the conducting states 
which produces an initial jump in the quantised conduc¬ 
tance from zero to 4e^//i. Our simulations provide no 
mechanism to explain the recovery of the first conduc¬ 
tance plateau as the confinement of the quantum wire is 
further reduced pAHIb] , 

The modulations in the transverse density calculated 
in this paper are open to experimental investigation. 
Using a weak magnetic field, electrons ejected from a 
quasi-one-dimensional wire can be focused into a quan¬ 
tum point contact (QPC). Variations in the transverse 
density, such as the double rows observed in Fig. 
would manifest themselves as variations in the current 
across the QPC as the magnetic field is varied [26] . 
Additionally, the tunneling of electrons from a quasi- 
one-dimensional wire into an adjacent two-dimensional 
electron system would contain signatures of the high- 
frequency transverse momentum components contained 
within the density modulations observed in this pa¬ 
per [23 EH] as well as other correlation and many-particle 
effects [29f[3T| . 

Devices of this kind have the potential for many appli¬ 
cations in quantum nano-electronics. Electrostatic and 
Coulombic interactions allow switching between one and 


two wires at length scales comparable to the surface gate 
features. The ability to dynamically define quantum 
wires at the micron scale could be an important compo¬ 
nent in integrated quantum architectures. For example, 
by overlaying the top gate over a small length of the split 
gate, a single wire could be split to form an Aharanov- 
Bohm ring with a tunable area, which could be used for 
high-precision magnetosensing. The spontaneous local¬ 
isation in the transverse direction demonstrated in this 
device could also be a precursor to Wigner crystallisation, 
which has applications in quantum information process¬ 
ing as a quantum memory or a spin-current mediator [32] . 


Additionally, the one-dimensional wires created by the 
double-row formation lie almost entirely underneath the 
split gates. In conventional one-dimensional quantum 
wires, electrostatic repulsion repels the electrons as far 
as possible from the split gate material. However, some¬ 
times it is desirable for the quantum channel to form 
closer to the surface gates. For example, proposals to cre¬ 
ate Majorana fermions in quantum wires [33U35] typically 
rely on the proximity effect, where Cooper pairs can tun¬ 
nel from a nearby superconductor into a one-dimensional 
wire which leads to bound electron pairs forming in the 
wire. The strength of this effect is exponentially depen¬ 
dent on the distance over which the Cooper pairs need 
to tunnel. The top-gate split-gate configuration studied 
here brings the electrons closer to the surface gate mate¬ 
rial and would enhance these effects, making it an ideal 
platform for studying proximity-effect induced physics in 
one-dimensional quantum wires. 
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In this appendix, we detail the composition of the semi¬ 
conductor heterostructures simulated in Sec. imi and El 
Eor both devices, the wave functions are fully enclosed 
within the simulation domain and we found no evidence 
of edge effects. 


1. Split-gate device 


Appendix A: Calculating the Dopant and Surface 
Charges 


The static charges consist of dopant and surface states 
which are thermally activated at high temperatures [36] . 
At 70K, these states are frozen in their equilibrium con¬ 
figuration and the fraction of occupied states is fixed. 
The resulting charge distribution provides the neutral¬ 
izing background for the electrons in the quasi-one- 
dimensional quantum wire and must be calculated before 
performing the main DET calculation. 

We determine this static charge distribution using a 
one-dimensional self-consistent calculation at T = 70K. 
The chemical potential of the GaAs surface is pinned 
in the middle of the band gap and we assume that it 
is also pinned at this value at the interface with the 
substrate. When the gates are initially grounded, the 
system is translationally invariant. Eor this simple one¬ 
dimensional calculation, Poisson’s equation is solved on 
a regular grid with a lattice spacing of Az = Inm. The 
energy of the dopant states is set to 5meV below the local 
conduction band minimum and the ionized dopant frac¬ 
tion is calculated using their thermal occupation. The 
density of electrons in the conductance band is calculated 
in the Thomas-Eermi semi-classical approximation. Con¬ 
vergence is achieved using the iterative method described 
in Sec. [TT] and the dopant density distribution at 70K is 
fixed. The carrier density at the GaAs/AlGaAs interface 
is then calculated at the base temperature using DET in 
the local-density approximation. The number of active 
dopants can be calibrated so that the two-dimensional 
carrier density is the same as that which is experimentally 
measured for a specific semiconductor heterostructure. 

As the semiconductor heterostructure is electrically 
neutral, the electric field above the surface of the sample 
vanishes. Therefore, the charge per unit area on the sur¬ 
face must bend the potential such that E = —eV(j) = ^ 
for z > 0. Erom Gauss’ law, this requires that 


<^surface — 



(Al) 


at T = 70K where crgurface is the charge per unit area of 
the surface. Laux et al. [2] indicate that there should be a 
20 nm depletion region of the surface charges around the 
metallic surface gates. We include this depletion region 
although its main effect on the quasi-one-dimensional 
quantum wire is to provide an offset in the gate voltages. 


The device simulated in Sec. mi is the same as the 
one found in the simulations of Laux et al. [2]. The 
device consists of two metallic surface gates at a volt¬ 
age Vg separated by a gap of 400nm on a semiconductor 
heterostructure with a 26nm GaAs cap above a 34nm n- 
doped layer of AIq. 26 Gao. 74 As with a donor concentration 
of Nd = 6.0 X 10^^cm“^. A lOnm AIq, , 26 Gao. 74 As spacer 
layer separates the dopants from a GaAs/AlGaAs inter¬ 
face against which the electrons are situated. This GaAs 
layer extends to a substrate which is 6 /im below the sur¬ 
face and contains a weak background p-doping with an 
acceptor density of Na = lO^^crn”^. Above the sam¬ 
ple is a 2 /im layer of vacuum with the natural boundary 
condition = 0 on the top surface. In these simu¬ 
lations, we set the density of charged surface states to 
^surface = —1-6 X 10^^cm“^ and we have included the 
20 nm depletion region of the surface charge adjacent to 
the gate metal, which gives a split gate voltage offset of 
-f-O.lV with respect to the original Laux et al. paper. 
The electron density at the GaAs/AlGaAs interface was 
calculated for T = 4.2K on a regular grid with Ay = 5nm 
and Az = Inm from pmin = — 200 nm to Pmax = 200 nm 
from the center of the split gate and = 90nm to 
^max = 64nm below the surface. 


2. Top-gate split-gate device 


The device modeled in Sec. IV is the one used by Ku¬ 
mar et al. in Ref. m- The two metallic surface gates are 
separated by a gap of SOOnm and the heterostructure con¬ 
sists of a lOnm GaAs cap above a 200 nm AIq.33 Gao. 67 As 
n-doped layer with a 75nm AIq.33 Gao. 67 As spacer. Below 
the spacer is the GaAs/AlGaAs interface against which 
the one-dimensional wire is defined. The GaAs extends 
to a substrate 2 /im below the surface and a 150nm layer 
of polymethylmethacrylate (PMMA) separates the sur¬ 
face from the top gate. The donor concentration in the 
AlGaAs dopant layer was = 4.3 x 10 ^^cm“^ which 
corresponds to a two-dimensional electron carrier con¬ 
centration of n2D = 1.0 X 10^^cm“^ at a base tempera¬ 
ture of T = 50mK. The electron density was calculated 
on a regular grid using the effective-band approximation 
described in App. The simulation regime for calcu¬ 
lating the electron density was from ^min = — lOOOnm 
to ^max = lOOOnm from the center of the split gate and 
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^min = —325nm, Zmax = “275nm below the surface. The 
lattice spacing of the grid was = 5nm and Az = Inm. 


Appendix C: The Effective-Band Approximation 

Solving the Kohn-Sham equation for a quasi-one- 
dimensional wire Eq. is computationally expensive, 
especially in the weak-confinement regime studied in 
Sec. |IV| where the wave functions can spread out in the 
^-direction. To reduce the complexity of the problem, we 
used the effective-band approximation of Stopa [37]. In 
this section, we show why this approximation works and 
the assumptions required to ensure its accuracy. 

Let us expand the wave function in terms of the eigen¬ 
states of the Hamiltonian in the z direction 

2^) = E (Cl) 

n,y 


where are the solutions to the one-dimensional 

Schrodinger equation 

+ =<^n{y)^Uz) (C2) 


for fixed y. As in Sec. |TT] the density dependence of the 
exchange-correlation potential is dealt with iteratively so 
we will assume that V{y,z) does not depend on n{y,z) 
for the purposes of the approximation. Inserting Eq. |C1| 
into Eq. j^for /c = 0, premultiplying by 
integrating over z gives 

\2^d^ ^ ) Cij,m{y) + fj,m{y) = Ejaj^m{y) 

(C3) 

where 


poo 

fj,m{y) = - ^ J CViz) 

/ oo 

C*(^) 

-c» 


^ daj^niy) 


E 


dy^ 


dz 


^ dy 

n,y 


dy 


-dz 


(C4) 


Taking ^ to be a parameter for the Hamiltonian in 
Eq. C 2 , Hellman-Eeynman theory [38] allows us to 
rewrite the integrals in Eq. |C4| as 


/. 

/ 


= IZo9yny,z)^V{zmz)dz 


— oo 
oo 


dy 


Smiy) - e„(y) 


(C5) 


c {dyViy,z)^Uz) + ‘^dyV{y,z)dy^y{z) - 2dyeniy)dy^Uz)) dz 




emiy) - en{y) 


for n m. We want to show that the term fj^m{y) vanishes for the electrons in layered, gate-defined semiconductor 
heterostructures. Eor Eq. |C5[ 


/ oo pOO 

dvV{y,z)ilp{z)il{z)dz< / \dyV{y,z)^^*{z)Cy{z)\dz 

-OO J —OO 

< \dyV{y,z)CK{z)f dz 

= \J{\dyV{y,z)\^'j^ 


(C7) 

(C 8 ) 

(C9) 


where the second inequality is due to the Cauchy-Schwartz inequality and (•)^ is the expectation value of an operator 
for a state in the eigenfunction We note that due to the orthogonality of the ^'^{z) 




(CIO) 


so using the Cauchy-Schwartz inequality on the various terms in Eq. |C 6 | gives 





d^^niz) 

dy'^ 


dz < 


/o 

\d^y{y^z)\ 



^{\dyV{y,z)\^)^ 

eniy)-ek{y) 

2 

dz 

em(y) - e„(y) 


(Cll) 


In semiconductor heterostructures, the electron wave functions are strongly confined in the growth direction so only the 
ground-state eigenfunction in the ^(-direction ^o(^) occupied and aj^n ~ 0 for n 7 ^ 0 and we only need fj,o{y) —> 0 . 




























11 


Any variations in the energy scales in the transverse direction along the heterostructure interface are orders of 
magnitude smaller than ei — eo- Additionally, for a typical heterostructure interface there are no degenerate states 
with \en{y) - ek{y)\ 0. Therefore, 


\eo{y) - en{y) \ » ^J(^\dyV{y,z)f^^, (^\d^V{y,z)\^'j^, dySniy) 


(C12) 


SO fj,o{y) -> 0 and Eq. C3 reduces to a one-dimensional 
Schrddinger equation with a band-dependent effective 
potential. The argument made here for the accuracy 
of the effective-band approximation extends to the two- 
dimensional electron systems in the original paper by 
Stopa. The density is given by 


This method significantly reduces the computational 
complexity of the electron density calculations which 
allows us to explore the low-density, weak-confinement 
regime. 


/ oo 

fiE,T)gME;sj) \aj,o{y)\^ \eoiz)fdE 

J 

(C13) 
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